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Abstract. We apply our recently developed theory of frequency- filtered and time- 
resolved TV-photon correlations pQ to study the two-photon spectra of a variety of 
systems of increasing complexity: single mode emitters with two limiting statistics (one 
harmonic oscillator or a two-level system) and the various combinations that arise from 
their coupling. We consider both the linear and nonlinear regimes under incoherent 
excitation. We find that even the simplest systems display a rich dynamics of emission, 
not accessible by simple single photon spectroscopy. In the strong coupling regime, 
novel two-photon emission processes involving virtual states are revealed. Furthermore, 
two general results are unraveled by two-photon correlations with narrow linewidth 
detectors: i) filtering induced bunching and ii) breakdown of the semi-classical theory. 
We show how to overcome this shortcoming in a fully-quantized picture. 
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1. Introduction 

At the heart of the modern theory of light lies the concept of the photon, the quantum 
of the electromagnetic field [2] . Regardless of the wave behaviour of light in a variety 
of contexts, ultimately, every measurement is accounted for by clicks on a detector. 
When the number of quanta is so large that the granularity of the field disappears due 
to coarse-graining, the limit of classical electrodynamics is recovered [3]. It is one of 
the wonders of quantum mechanics that the light field still exhibits wave-like properties 
when tracked at the single-particle level. Most popularly, the interference patterns in a 
double slit experiment is reconstructed by the accumulation of a large number of single 
clicks, each observed in isolation [4]. 

Usually, even in the regime of a small number of quanta in the field, the experimental 
observables are obtained by averaging over a large enough repetition of the experiment. 
For instance, a luminescence spectrum is obtained by integrating the signal over long 
times, such as to obtain a continuous distribution that describes the probability of 
emission of each photon at a given energy. An example for such a one-photon spectrum 
(IPS) S(lj) is given in Fig. [l](a) for the simplest conceivable case, that of a single mode 
emitter. The lineshape is simply a Lorentzian with a linewidth that is given by the 
inverse lifetime, corresponding, in the time domain, to an exponential decay of the 
excitation: the emission is memory less and has the same probability to occur at any 
given time, as befits a quantum particle. This results in a spread of energy of the emitted 
photon. Arguably, however, this lineshape is also typical of a classical system, in the 
form of a damped resonance. 

To access the dynamics of photons at a quantum level, with no counterpart from 
a classical theory involving a continuous field, it is necessary to invoke A-photon 
correlations. An average is still taken over a large number of single-detection events, 
but the latter now involves N correlated clicks in each shot. In this text, we will 
focus on the simplest case where N = 2. There is an ample literature on two-photon 
correlations [5j. Typically, photons are detected as a function of time (and/or other 
commuting variables like space) in photon counting configurations such as the celebrated 
Hanburry Brown-Twiss setup [6]. The specificity of our approach is that we study 
correlations in the frequency domain, i.e., we compute correlations between a photon 
detected at an energy uj\ and another at energy U2. The conjugate time information 
is still present through the linewidth of the frequency windows. Real time can also 
be included, computing correlations between a photon of frequency uj\ detected at 
time t\ and another of frequency uj 2 detected at time t 2 ? although we will limit here 
to simultaneous detection for brevity. This is a generalisation of the conventional 
(color-blind) theory. The two-photon spectroscopy of a quantum emitter is obtained 
when such correlations are computed (theoretically) or measured (experimentally) for 
all the energies spanning the range of emission of the system. Experimentally, the 
implementation is straightforward, as it essentially consists in inserting filters between 
the system and the conventional photon-counting setup. The technique is increasingly 
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popular in the laboratory [Zl El El HDl CD] although its interpretation relies chiefly on 
physical grounds or numerical simulations. Theoretically, almost all the effort has 
been dedicated to resonance fluorescence [T21 [13l El HSl HE], owing to the difficulty 
of actual computation for other systems, prior to our input pQ that yields exact results 
for arbitrary systems. 

In this text, we present a comprehensive analysis of a variety of fundamental 
quantum optical systems and show the extremely rich dynamics and the novel insights 
that such frequency-resolved two-photon correlations unravel. In Section [21 we give a 

(2) 

brief overview of our theory to compute two-photon spectra (2PS), denoted g], (^1,^2). 
The notation refers to the Glauber precept that places photon correlations as the pillar 
of quantum optics [5]. The mode which photons are collected from will be appended 
in square brackets for clarity but only when necessary to avoid confusion. The full 
theory for A^-photon correlations and/or arbitrary time delays is laid down in Ref. pQ. 
In Section |3j we address the simplest case of a single mode: an harmonic oscillator or a 
two-level system. These serve as the first illustration of our theme that the 2PS reveals 
a considerably richer picture of the dynamics of emission for an otherwise identical 
IPS lineshape, depending on both the mode statistics or the nature of its broadening 
(homogeneous or inhomogeneous) . In Section [ij we analyse the more complex dynamics 
of coupled modes. Here again, even though the low-pumping regime results in identical 
IPS lineshapes, the coupling of two harmonic oscillators, two two-level systems, or one 
harmonic oscillator to a two-level system (Jaynes-Cummings model), show distinct and 
peculiar features in the 2PS, unravelling the different nature of the three couplings. 
In the latter cases, of great importance in cavity-QED, we thereby identify an entire 
class of two-photon emission processes, to which one is completely oblivious through 
one-photon spectroscopy. In Section [5| we address the high excitation regime of a 
two-level system, considering two ways of driving it coherently: with the standard semi- 
classical model [17J, that is, through the coupling to a classical term in the Hamiltonian, 
describing an ideal laser, and with a fully-quantized model PS] where a very large 
number of photons builds up a coherent field in a cavity. Both provide the expected 
Mollow triplet lineshapes in the IPS but only the full quantum model provides a physical 
2PS. In Section [6| we conclude and discuss further directions. 

2. Theory of one- and two-photon spectra 

We have recently shown how to compute efficiently and easily two-photon correlations 
between photons of frequencies u\ and cj 2 through the coupling of two "sensors" to the 
open quantum system of interest [lj. In the limit of vanishing coupling of the sensors to 
the mode which photons correlations are to be measured, the sensors intensity-intensity 
cross-correlations recover exactly the frequency-resolved two-photon correlations defined 
in the literature through integrals that are, however, too cumbersome to compute 
directly in most cases [19]. In the case of zero delay between the photons, the normalised 
correlation function, or two-photon spectrum (2PS), is therefore obtained simply in terms 
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of the population operators of the sensors ri\ and 712, as: 

(2) / \ r jnin 2 ) , 
grur ^,u; 2 ) = hm^— y (1) 

where e represents the coupling to the sensors, small enough not to disturb the system, 
and I\, T 2 are the sensors decay rates that provide the detectors linewidths or, in other 
terms, the respective widths of the frequency windows. The case of equal frequencies 
g^\cj,cj) measures the photon-statistics of photons detected in a frequency window of 
Lorentzian shape of linewidth Y. This corresponds to applying a single filter or to the 
effect of the detectors resolution in a colour-blind measurement, and is thus also of prime 
interest. 

In practical terms, the calculation can be performed in two alternative ways, as 
discussed in Ref. pQ. The first method relies in the actual inclusion of the sensors in 
the dynamics of the system, supplementing the Hamiltonian and the master equation 
with the corresponding terms. The main drawback is that the size of the Hilbert space 
is increased by a factor 2 2 (and the density matrix by a factor 4 2 ). However, the 
computational operations are simple. In one of the possible computational approaches 
to the problem, one merely needs to solve the steady state of the full system (a set of 
homogeneous linear equations), which only requires the inversion of a matrix [20] . An 
alternative method, presented in detail in the supplemental material of Ref. [lj, consists 
in using the formal semi-analytical expressions for the populations and cross correlations 
between the sensors to leading order in the couplings, e 2 and s 4 respectively. Then, the 
relevant quantities (ni), (n 2 ) and (nin 2 ) can be obtained through these expressions (that 
are provided below) and the master equation of the system only with no increase in the 
numerical complexity. In order to implement this method, one first needs to determine a 
vector of the system operators whose average are needed to compute the coupling to the 
sensors. Calling O the annihilation operator of the mode of interest, which correlations 
one wishes to compute, such a vector can be written in the form v = (1, O, 0\ 0^0 . . .) 
(with subsequent terms depending on the systems and its dynamics). The observables 
of interest for the sensors cross-correlations are specified by the mean values of this 
vector: (v) = (1, (O), (O^), (0^0), . . .). A regression matrix M can be obtained from 
the master equation of the system [20j which rules the dynamics of these correlators 
according to d t (v) = M(v). The steady state (if it exists) is obtained in the limit of 
infinite times: (v ss ) = lim^oo e Mt (v(0)), for any initial state v(0). Next, one builds 
two re-ordering matrices, T±, which, when acting on the vector (v), introduce in all 
correlators an extra operator for T + and an extra operator O for T_, keeping 
normal order in each case. That is, T + (v) = ((O f ), (0+0), ((O f ) 2 ), ((O f ) 2 0), . . .) and 
T_(v) = ((O), (O 2 ), (0^0), (O^O 2 ), . . .), respectively. With these matrices, one obtains 
the populations of the sensors to leading order in e, as the first element [. . .]i of the 
vector: 
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Regarded as a function of the frequency u>j, this produces the IPS as measured by 
a detector of linewidth Tj. Similarly, the cross-correlations that provide the 2PS is 
given by: 

2e 4 r„ -1 

(™1™2) = ~ & 
J- 1 + J- 2 



M + (-iw2-ri-^)i 

~ L X T, -— ^— T_ + T_ _1 



M-raV + M + (-^!-^)l M + (iwi-^)l 

1 't_ „ ~ l fr— — T + + T H - J 



M + (^ 1 -iw 2 -^±^)lV M + (M-^)1 + + M + (-iw 2 -^)l 

zl T_ f — - I — - ^T_lv ss 

+ M + (-icui -icu 2 - \M + (-«j 2 -^)1 M + (-iwi-^)l/ / 

+[1 o 2] , (3) 

where [1 o 2] means the interchange of sensors 1 and 2, that is, permuting u\ o cj 2 
and Ti T 2 everywhere. 

The advantage of this second method is twofold. First, it is very useful when the 
Hilbert space is small as it can lead to closed-form analytical expressions, as we show 
in the next Section. Second, it may also be numerically advantageous as the matrices 
involved correspond to the original (smaller) Hilbert space of the bare system, in the 
absence of sensors. The computational price to pay is in a higher number of matrix 
operations (eleven different matrix inversions and numerous multiplications for (774712)). 
In some situations, the advantage of a smaller Hilbert space dominates. In others, it 
is more convenient and straightforward to explicitly include the sensors and solve for 
the enlarged system, with a similar overall numerical efficiency. We used both methods 
indistinctively to obtain the results presented in the rest of this text. 

From now on, the linewidths of the detectors will be taken equal, Y\ = T 2 = T, 

( 2) 

for the sake of simplicity. We will thus simply write g-p (^1, ^2)- The width T of 
the frequency windows is a fundamental parameter that cannot be dispensed with, 
unlike the IPS case where an ideal detector is usually assumed in theoretical works. 
In the 2PS, T provides the uncertainty in the frequency (r) and time of detection 
(1/r) [21], as required by the Heisenberg principle. In the limit of very broad filters 
(r 00) the standard (color-blind) second order coherence function at zero delay is 
recovered: limr^oo 9y > (^i? ^2) = 9^ • The opposite limit of very narrow filters (r — » 0), 
results in a systematic tendency regardless of the underlying system, namely photons 
are uncorrelated, and those with the same frequency are thermally bunched. This can 
be understood on physical grounds: in order for the detectors to provide high precision 
in frequency (r —} 0), their interaction time with the system has to be long (1/r 00), 
so that the collected photons correspond to all the possible times in the dynamics. 
This leads to an apparent uncorrelated statistics: lim r ^ dr^i^u ^2) — 1 provided that 
uji ^ UJ2- The limiting value becomes lim r ^ u ) = 2 in the case of detection 

at equal frequencies as there are 2! ways that two indistinguishable photons can be 
collected by two different detectors. We refer to this correlation tendency to bunch 
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Figure 1. (a) The one-photon spectrum (IPS) common to a single harmonic 
oscillator and a two- level system, and the corresponding two-photon spectra (2PS), 

(2) 

#P [O](o;i, CJ2), for (b) an harmonic oscillator and (c)-(d) a two- level system, in 
units of the total broadening To- In (d), the emitter is subject to pure dephasing, 
70 = 0.75r a , with T a such that the IPS remains the one shown in (a). Beyond the 
bunching on the diagonal, common to all 2PS, clear deviations are observed depending 
on the type of emitter: the harmonic oscillator is structureless while the two-level 
system concentrates its antibunching in a butterfly-shaped region of non-overlapping 
frequencies. Dephasing favours antibunching. The color code is: blue; 1 white; 2 
red. 



when uj\ ~ cj 2 and T is small (as compared to the relevant linewidths in the system), 
as indistinguishability bunching. This effect has been observed experimentally in the 
filtering of a single mode laser [22] . 

3. Single mode emitters 

Two-photon frequency-resolved spectroscopy is such a nascent field of optical 
characterization that even the trivial systems need to be investigated. Namely, our 
starting point is the free mode (or, in the quantum optics terminology, "single mode"), 
which Hamiltonian simply reads: 

H = ujqO^O . (4) 

There are two fundamental possibilities for the operator O, namely, it can be the 
annihilation operator of an harmonic oscillator, in which case we denote it O — a, 
or it can be that of a two-level system, in which case we denote it O — a. Their 
respective quantum algebra is given by the bosonic commutation rule {ao) — a^a = 1) 
and the fermionic anticommutation rule (aa^ + aV = 1). The dynamics of the density 
matrix p of the emitter — including decay (required to bridge with the external world 
where the measurement is performed) and a continuous incoherent pump (to populate 
the system) — is given by the master equation 

d tP = i[p, H ] + ^£o(p) + ^y£ot(p) , (5) 
with the Lindblad super-operator Co{p) = (20 pO^ — O^Op — pO^O). 
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The IPS, provided by Eq. @, is the same for both systems — a, a. It is plotted 
in Fig. [jja). It is a Lorentzian lineshape with linewidth Tq (inverse mode lifetime) 
given by T a = j a — P a for the harmonic oscillator and T a = j a + P a for the two-level 
system [23J. 

The two-photon correlations, however, are dramatically different, since the 
harmonic oscillator displays bunching, g^[a] = 2, while the two-level system displays 
the exact opposite behaviour, antibunching, g^[cr] = (the mode which statistics is 
measured is here denoted explicitly in square brackets to lift the ambiguity). The 2PS 
also exhibits structures in the frequency correlations that are qualitatively different. For 
such simple cases, Eq. ^ can be solved exactly and analytical formulas for the 2PS 
obtained: 

with a common expression regardless of the mode O = a, a spelt out above, and a term 
(?p 2 ^[0] specific to each case, given by: 

#V] = 2r(i + £) , (7) 

for the harmonic oscillator (note that it is frequency independent), and by 

(3r+rcr 3r+rcr \ 

for the two-level system. The single filter case g\?\u,oj) shows that for the two- 
level system, correlations are maximum when filtering the peak itself, uj = 0, and 
that a large overlap of the peak is needed to recover a good antibunching of the 
two-level system (the ideal result is recovered in the limit of broadband detection: 
limr^oo <7r^ M O^i > ^2) — fl^M — 0). The photon-counting correlation of the filtered 
peak indeed reads ^p 2 ^[a](0, 0) = 2(T a /T)/{3 + T^/r), e.g., the filter linewidth must be 
over sixty-six times that of the emitter to reach the 1% level of accuracy, an awkward 
requirement in practice. For this reason, post-processing of the coincidences through 
deconvolution of the spectral response is important to characterise fairly a quantum 
emitter. 

For the harmonic oscillator, on the other hand, g^[a](uj,uj) = 2 regardless of the 
filtering window T and of which part of the peak is detected (maximum or any point in 
the tail). This is a manifestation of the classical character of the harmonic oscillator: it 
has no local information, any part behaves like the whole and one cannot tell apart from 
a filtered window whether the information is of a microscopic or macroscopic nature. 
In contrast, the two-level system has an energy scale and the information is localised in 
the energy window proper to the dynamics of the emitter. 

The full 2PS of these systems further reveals such fundamental aspects of these 
emitters. In Fig. [TJ we plot Eq. ^ in the case where the detector linewidth matches 
that of the IPS peak, V = Tq- In this and all subsequent density plots of the 2PS, we 
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use a color code where red corresponds to bunching (g^ > 2), blue to antibunching 
(g^ « 0) and white to non-correlated emission (g^ = 1). We are more concerned 
in this text with the qualitative patterns that emerge in two-photon spectroscopy than 
quantitative results (the optimisation of given correlations is straightforward but would 
bring us to a too voluminous discussion) . 

The common characteristic between the harmonic oscillator and the two-level 
system is the indistinguishability bunching line, for frequencies that are indeed 
indistinguishable within the detector line width \uj\ — c«J2 | < T. This accounts for the 
diagonal feature on all plots. Apart from that, the harmonic oscillator correlations 
lack any structure and are always above one, corresponding to the expected bunching, 
whereas the two-level system correlations mainly assume values below one, the expected 
antibunching, furthermore in a nontrivial configuration of detection C with different 
frequencies that attempt to maximise the overlap with the IPS without entering the 
indistinguishability bunching region: 



The largest antibunching is thus obtained on the antidiagonal. In practical terms, for a 
given detector linewidth, it is therefore better to perform photon coincidences between 
the left and right elbows of the IPS, with a slight shift of both windows away from the 
center, rather than between photons both coming from the central peak. One can easily 
compute where exactly to place the filters by evaluating g^a]^, — uS) with Eqs. (|6j) 
and M. For instance, in the case Y = T a , the minimum goes down to 2(\/2— 1)/5 ~ 0.16, 



for uji = —uj 2 = y (15\/2 — 4)/31 ~ 0.75 IV This is the closest that the frequencies 
can get to the central one without suffering from the indistinguishability bunching. It 
is a considerable improvement, at the small cost of a reduced signal, on measuring the 
correlations from the central peak with a detector of the same linewidth, that provides 
an antibunching of 0.5 only. When both detection windows are far from the peak, the 
emitter looses its quantum character and behaves like its classical counterpart, exhibiting 
bunching of its photons on the diagonal and uncorrelated emission otherwise. 

We complete the study of the correlations of a single two-level system emitter by 
characterising the effect of inhomogeneous broadening. This is relevant, for instance, 
in semiconductor systems (quantum dots) where the solid state environment induces 
fluctuations and pure dephasing on the levels. In this context, filtered two-photon 
correlations have already shed light on the timescales and origin of fluctuations in 
quantum dots [10J. Theoretically, inhomogeneous broadening is introduced in the 
dynamics through another Lindblad term of the form -y/V<r(p) EH EH]. Pure dephasing 
occurs at the rate 7^, diminishing coherence in the system without affecting directly the 
populations. Its effect on the IPS is merely to increase the linewidth of the peak, 
T a = 7(7 + ^ + 7^- For this reason, it is typically difficult to measure a radiative lifetime 
from spectroscopy and time-resolved measurements are usually invoked for that purpose. 
For a given IPS, however, the 2PS changes quantitatively, which allows radiative and 
non-radiative contributions to be measured directly from spectroscopy measurements. 



C = 




(9) 
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The 2PS spectra without (c), and with (d), pure dephasing are compared in Fig. [TJ 
While the statistics of the two-level system as a whole is not affected by 7^, its 2PS 
structure is. Namely, pure dephasing extends the condition of antibunching Q to a 
wider range of frequencies and enhances the ant icorrelat ions. This counter-intuitive 
result is a manifestation in 2PS of the profitable effect of dephasing for some quantum 
correlations [26l [27] . 

4. Coupled single mode emitters 

The bare emitters already display a rich structure of two-photon correlations. We 
now show that the next simple system, that of two emitters coupled linearly, also 
exhibits two-photon correlations of unsuspected complexity. Two coupled single- 
mode emitters account for, or are at the heart of, a large gamut of quantum optical 
systems. We will consider three paradigmatic cases, subject of intense theoretical and 
experimental studies: two coupled harmonic oscillators [28l[29], two coupled two-level 
systems [30], EU E21 [33] and the Jaynes-Cummings model [3H |35j EJ EHl EH [38] where 
one harmonic oscillator is coupled to a two-level system. All of them are described by 
the following Hamiltonian: 

H c = u 0l O\Oi + u 02 0\0 2 + g(0\0 2 + 0\O{) , (10) 

where, as in the previous section, Oi^ describes two modes that can be either an 
harmonic oscillator or a two-level system. The dynamics of the system is completed 
with the inclusion of decay and pumping through the master equation: 

d tP = i[ P , H C ] + [^fc Ql + lf£o 2 + ^fC o{ + ^C ol ]( P ) . (11) 

In the rest of the text, for the sake of brevity, we will limit to the situation where the 
first mode 0\ is detected and the other mode O2 is pumped (Po 1 = 0). This reduces the 
number of parameters and remains close to the experimental situation in cavity-QED 
where typically the emitter is incoherently pumped and the system is observed through 
the cavity emission. 

In Fig. [2j we show the 2PS of the Oi-mode for different values of r )o 1 , from Ag (weak 
coupling) to O.lg (strong coupling), fixing the lifetime of the second mode to a small value 
7o 2 = O.OOlg, also on the basis of typical experimental values. The detectors linewidths 
are set equal to that of the emitting mode, T = 70^ as previously. In the linear regime, 
when Pq 2 —> 0, all the models converge to the same IPS (first row), whereas their 
2PS exhibit very different structures depending on the nature of the emitters. When 
the losses overcome the coupling, i.e., in the weak-coupling regime, the IPS converges 
to that of the bare mode 0\. In the case of the coupled harmonic oscillators and the 
coupled two-level systems, the 2PS has the form of the single mode, cf. Fig.[l](a) and (b), 
respectively. However, in the case of the Jaynes-Cummings model, although the cavity 
mode is detected, a statistics reminiscent of the two-level system is observed, even in 
weak-coupling. This hints at the value of such systems for quantum applications, such 
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Figure 2. First row: IPS that are, in the linear regime, identical for coupled 
harmonic oscillators, coupled two-level systems and the Jaynes-Cummings model. 
Second/third/fourth rows: corresponding 2PS, that are, on the contrary to the IPS, 
markedly different. In the coupled two- level systems and the Jaynes-Cummings model, 
a logarithmic scale has been used to improve the visibility of all the features. Columns 
(a)-(d) correspond to increasing the quality of the coupling as r )o 1 /g = 4, 2, 1, 0.5, 0.1. 
In all cases, T = O.lg, Po x = and 7o 2 = O.OOlg. 

as a single-photon sources: the (classical) cavity enhances the emission and collects 
the light while still retaining the valuable quantum features of the (quantum) emitter. 
Note that while juxtaposing the two filters on both sides of the central peak was only 
an improvement in the case of the bare two-level system, it is here mandatory, as the 
central peak is strongly bunched. 

As the quality of the coupling increases (7^ — >> O.lg), the system enters in the 
strong-coupling with new dressed states, or polaritons, emerging. We call p± the 
annihilation operators for the upper (+) and lower polariton (— ). The IPS, still identical 
for the three coupling models, turns into the so-called "Rabi doublet", with splitting 
2R = 2a/ g 2 — [(701 7o 2 )/4] 2 between the polaritons (the splitting observed in PL is 
more complicated and is given in Ref. [39]). The 2PS, however, present, again, different 
features that we discuss in turns. 
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4-1. Coupled harmonic oscillators 

The coupled harmonic oscillators, that we denote a and 6, are characterised by 
bunched correlations between any two frequencies, as expected due to the bosonic 
character of both modes. Similarly to the case of a single harmonic oscillator, the 
color-blind statistics always remains thermal, for both the bare and dressed modes: 
g^[a] = g^[b] = g^[p + ] = g^[pJ\ = 2. The cross correlations, however, depend on 
the system parameters: 

(2) (qtq&tfc) _ Pg/(n a ) + P b /(n b ) 

{n a ){n b ) r a + r 6 

(at resonance) and is always between 1 and 2. In the strong-coupling regime, g^[a; b] 
1, i.e., the two modes become uncorrelated. The correlations between polaritons, 
#( 2 )[p + ;]9_], are also close to 1 as long as polaritons are well defined (well separated 
spectrally), due to the formation of their independent polaritonic dynamics. 

Coming back to frequency-resolved correlations, the 2PS of the mode a (2nd row 
in Fig. [2]) presents, in strong-coupling, two hyperbolic dips where g^\^i^2) drops to 
its minimum value of 1. The hyperbolas, given by the equation 

uj r uj2 = -g 2 , (13) 

corresponds to the detected frequencies (^1,^2) hitting the polariton branches which, 
as a function of detuning A, read lj±(A) — A/2 ± ^ g 2 + (A/2) 2 . When one detects a 
photon at uj\ , the classical system does not provide information on whether this photon 
comes from the main peak itself (where the dressed state is located and the mode emits 
predominantly) or from any other point on its tail. This photon is detected as if emitted 
by the polariton being at precisely this energy. It has thus the fitting correlation of the 
coupled oscillators with another photon detected at the energy of the other polariton 
with the same detuning as the would-be polariton accounting for the first photon. The 
2PS thus evidence the existence of the full polariton dispersion via a change in the 
statistics, even though the modes are at resonance in the IPS. This shows again the 
richer discriminative power of two-photon spectroscopy. The 2PS anticrossing is wider 
than the IPS one given that the distance between vertices (±g, =FflO is 2^/2$ instead of 
2g. The renormalisation by \[2 is due to the fact that the information is carried at the 
two-photon level instead of one in the photoluminescence. 

When the modes are detuned, the polariton correlation interferences still form an 
hyperbola, 

(u>! - cj hyp ) (u 2 - cj hyp ) = -y h 2 yp , (14) 

but its center, c^hyp, and vertex, Vh yp , depend on the set of relevant parameters of 
the system, namely {P a) P b) j a , 75, uo a) uj b ] ) as shown in Fig. |3| In order to gain 
further insight into this dependence, let us assume that polaritons are well defined and, 
therefore, the IPS is essentially the sum of two Lorentzian peaks associated with each 
polariton, S[a](uj) w L a _(uj) + L°^{uj). We can neglect the dispersive contributions to the 
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Figure 3. 2PS for coupled and detuned harmonic oscillators. The bare states are 
at frequencies uo a = and = 5g (dashed lines). The total decay in the system is 
kept fixed to j a + 75 = O.lg while j a and 75 vary as indicated above each plot. Other 
parameters: P a = 0, Pb <C 7 a ,6 and T = O.lg. 



spectrum [40] thanks to the negligible overlap between these peaks and assume complete 
uncorrelation between the polaritons, g^[p + ]p^] « 1. In a single particle picture, one 
can define a frequency dependent operator for mode a, a{u) = y 'Ll {uj)p- + y / L^(cJ)p + , 
using the relative spectral weights, 

11 7±/2 



(15) 



7T ( 7± /2)2 + ( W - W± ) 2 ' 

that take into account the quantum state of the system through the coefficients l± [40j . 
The parameters 7± and u± are the full-width half-maximum and positions of the 
polariton modes, p±. The IPS is well approximated by (a^(uj)a(oj)). In a similar 
way, (a^ (ui)a^ (oj2)a(oj2)a(oJi)) , captures some of the features of the exact results in 
Fig. [3| The correlations computed in this way are less well reproduced than the IPS as 
this approach neglects the multi-photon dynamics. However, due to the linear nature 
of the system, they produce interferences between the well defined polaritons with the 



same hyperbola, Eq. (14), specified by: 
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In Fig. [3j we exemplify how the quantum state affects the anticrossing of the 2PS 
for a detuned situation, where cj_ « u) a = and uu + « uj^ = 5g. For simplicity, only 
one mode is pumped in the linear regime (P a = 0, P5 <C 75). The total decay of the 
system is fixed (j a + 75 = O.lg) whereas the relative value of the two rates is varied 
to alter the quantum state of the system. In the asymmetric case where 75 = 0, Fig. 
3](a), the normalized weight is 1+ ~ and, in accordance with Eq. (16), the center of 
the anticrossing ~ u a with a vertex Vh yp ~ g. In the symmetric situation where 
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7 a = 75, Fig. [3](b), both weights are equal, 1+ w P , yielding a balanced position of the 
anticrossing between the two modes, cJhy P ~ ; and a larger splitting given by the 

same Rabi splitting as the IPS, Vhyp ~ \J 9 2 + {w~a ^b) 2 : /4. Finally, in Fig. 3](c) where 
7 a = 0, the quantum state with P w 0, yields cjhyp ~ cj^ and the same vertex as in the 
case of Fig. ga), Vhyp « 



Coupled two-level systems 

The coupled two-level systems (3rd row in Fig. [2]) differs from the coupled harmonic 
oscillators in much the same respect than the two-level system does from the harmonic 
oscillator. There are now, unsurprisingly, regions of antibunching and a richer pattern 
of correlations. 

As the system enters strong-coupling, a grid of two vertical and two horizontal 
lines of antibunching appear, corresponding to one frequency fixed at one of the peaks 
at ±R and the other frequency scanning through the IPS. The stronger antibunching is 
obtained at their crossing. There is also the diagonal of indistinguishability bunching. 
Finally, a neat feature revealed by the 2PS is the strong superpoissonian antidiagonal, 
when uji + UJ2 = 0. This is the first occurrence of an important and recurrent 
pattern in two-photon spectroscopy: these correlations correspond to direct two-photon 
relaxation from the doubly occupied state =4 |0,0), through a virtual (off- 

resonant) intermediate state. As the energy of this auxiliary step is not constrained, 
the feature appears on an entire line. Only thanks to the coupling can one two-level 
system effectively emit two excitations at the same time (within the detector time 
window 1/r). We use the denomination of "leapfrog" to label this type of relaxation, 
as the intermediate rung is not populated and the system effectively jumps over it 
by emitting two photons simultaneously. It becomes more prominent as the detector 
width, r, decreases and the time uncertainty of the measurement increases. This is a 
fundamental process that is common to all quantum nonlinear systems. It supports 
the idea that the emission of two identical photons can be enhanced (and selected) by 
placing such systems in a high-Q cavity j4TJ |42j [43]. In the simple system that is the 
two-coupled two-level systems, it appears in isolation. We discuss it in more extent in 
the next paragraph where it appears along with similar processes, making sharper its 
own features. 



4-3. Jaynes-Cummings model 

In mixing the types of emitters, one harmonic oscillator (cavity) and one two-level system 
(atom, quantum dots, superconducting qubit. . . ), a much richer dynamics is unravelled 
by the 2PS (4th row in Fig. [2]). In strong-coupling, a ladder-type level structure (shown 
in Fig. i]) arises thanks to the nonlinear splitting 2y/ (\fng) 2 — [(j a — 7 cr )/4] 2 between 
the states with n excitations [20j . The one-photon relaxations between these states give 
rise to a rich multiplet spectrum [Sj, and to a considerably richer set of two-photon 
relaxations, that we now discuss in details. 



Two-photon spectra of quantum emitters 



14 



3 u 



-2 



















(a) 


(b) 

\ N \V/ / s 

\ V\ X v* ' ' 



















































































































































-2 







-4 








Figure 4. Patterns of two-photon correlations in the Jaynes-Cummings model up 
to the third rung of excitation (low pumping regime), (a) The grid of horizontal 
and vertical lines (solid black) corresponds to correlations between real states, when 
one frequency is pinned at a transition in the ladder that is correlated with emission 
anywhere else in the system resulting in resonances when matching with a cascade 
(bunching) or an incompatible process like emission from the other type of polariton 
(antibunching). (b) The grid of antidiagonal lines (dashed green) corresponds to 
leapfrog processes where two-photon transitions take place from one rung of the ladder 
to another one two steps below, jumping over the intermediate rung. The grid of 
diagonal lines (except the central one, dashed-dot blue) correspond to polariton-to- 
virtual-state antibunching. The "indistinguishability bunching" (diagonal red line) is 
common to all 2PS. These two panels overlapping emerge clearly in the best systems, 
cf. Fig. |2jc)— (e) , up to the second rung, and Fig. [5jc),(d). Panels i-iv on the right are 
sketches of the underlying processes responsible for these correlations (with the same 
color code). 



The simplest and most straightforward type of two-photon correlations is the one 
that precisely follows from sequences of two one-photon relaxations, with each photon 
emitted from one state of the Jaynes-Cummings ladder to another. Depending on the 
type of relaxation, a given type of correlation is expected. For instance, at low pumping, 
detecting emission from the two types of polaritons in any given rung of the ladder leads 
to antibunching, as the de-excitation went one way or the other. Similarly, detecting 
emission in cascades from the ladder in consecutive events (as depicted in Fig. leads 
to bunching. Such processes are the extensions to all possible rungs of the ladder of the 
correlations observed with the coupled two-level systems. In fact, when the system is 
not in very strong-coupling, and thus makes it difficult to resolve the dynamics of the 
higher rungs, the 2PS of the Jaynes-Cummings resembles that of the coupled two-level 
systems, with four antibunching dips appearing at uo\ or cj 2 = ±i?, corresponding to the 
situation depicted in Fig. [4jL 

As the system gets deeper into strong-coupling, more such correlations from 
transitions between different rungs are revealed. They are given by u\ = E n ±E n _ 1 and 
uj 2 free (vertical lines) and their counterpart uj\ <r+ uj 2 (horizontal lines). All together 
they constitute the grid of horizontal and vertical lines shown in black in Fig. [4](a). 
This pattern is neatly reproduced by the actual 2PS (see Fig. [5| for the best system) 
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Figure 5. IPS (first row) and 2PS (second row) for Jaynes-Cummings, going from 
vanishing pumping in (a), (c), to the nonlinear regime in (b), (d), P a = 0.05g. Second 
rung transitions appear in the IPS inside the Rabi doublet and new features involving 
the third rung of excitation appear in the 2PS. 

with enhancement/suppression of the correlations at every crossing point depending 
on whether the combination of frequencies corresponds to consecutive/non-consecutive 
processes pp. 

There is also the now familiar diagonal of indistinguishability bunching, in red 
in Fig. |i](b), and other diagonal and antidiagonal features, better resolved in the 
best systems. These correspond to emission through virtual states during the short 
interaction time with the detector, ~ 1/r, and have, to the best of our knowledge, not 
been reported up to now. The two antidiagonal lines at uoi + u) 2 = 2uo a ± E 2 ~ ±V%g 
(dashed green lines in Fig.[i](b)) are the leapfrog processes already encountered with the 
coupled two-level systems: a two-photon transition from the second rung directly to the 
ground state, transiting through a virtual intermediate state which unspecified energy 
allows this process to be detected at any point on a line. This process is depicted in 
Fig. [pn. The correlation for such a direct two-photon emission is, naturally, bunched. 

Increasing the excitation power, the higher rungs get populated and the additional 
leapfrog process |3±) =4 |1±) can be observed. The four satellite antidiagonal lines that 
result, satisfying u)\ + uj 2 = 2uj a ± (E 3 — E\) or ±(£3 + £1), are plotted in Fig. |i](b) and 
materialised in the 2PS in Fig. [5](d). These lines are more visible in some regions (of 
strong bunching) but upon closer inspection, a deviation of the statistics due to such 
processes is realized everywhere. 

A feature not observed in the previous cases is the secondary diagonal lines at 
UJ1—UJ2 = ±2i?, that meets with the antibunching dips at the first rung resonances. They 
are plotted as dotted-dashed blue lines in Fig.[i](b). In the actual 2PS, these are globally 
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weaker and more clearly distinguished in some regions only, like the leapfrog lines from 
the higher rungs, and are also realized everywhere if looked at closely. Unlike leapfrog 
processes or indistinguishability bunching, however, they do exhibit an antibunching 
tendency as compared to the surrounding region (even if the overall correlation remains 
bunched). They correspond to the anticorrelated emission from both polaritons of the 
first rung into the same virtual final state, as depicted in Fig. ^v. These lines become 
more pronounced when the mode under observation is also the one that is excited (as 
stated previously, we consider the situation where one mode is excited and the other 
one is observed). In such a pumping configuration, strong polariton-to- virtual-state 
correlations also appear in the coupled two-level systems (not shown) while, as seen in 
Fig. [2j third row, they are absent in the other pumping-detection scheme. This suggests 
that the interference between relaxations from |+) and |— ) plays a role. They are thus 
not peculiar to the Jaynes-Cummings model but to any strongly coupled, nonlinear 
system. Higher excitation, makes higher order diagonal lines appear, corresponding to 
higher rung dressed states, such as u)i — u) 2 = ±2i? n . Such diagonals from the second 
rungs are also visible in Fig. [5](d). In this figure, the effect of increasing pumping is 
shown to "melt" the pattern of correlations as more real states are populated and more 
features distributed over the 2PS. 

All these results show that two-photon spectroscopy unfolds the mechanisms of 
relaxation that are not apparent in photoluminescence. This can be used to characterize 
a system or exploit quantum correlations. 

5. Lasing and the Mollow triplet 

We have just discussed the effect of increasing the excitation power in the case of 
the Jaynes-Cummings model: higher rungs get populated and give rise to additional 
relaxation processes, either from real states that involve photon by photon de-excitations 
or through two-photon emission that involve intermediate virtual states. Further 
increasing P a may bring the Jaynes-Cummings system into the lasing regime [151 1^6] 
where the optical field becomes coherent and g^[a] — >■ 1. The cavity IPS reduces 
to a single line that narrows with increasing occupancy of the mode as 7l ~ 
g 2 / (2^ a (n a ) 2 ) [47J, where (n a ) ~ P a /(2~f a ). The inverse quantity, 1/7l, corresponds 
to the coherence time of the laser. More interesting is the IPS of the two-level system 
(or quantum emitter), which converges to a Mollow triplet [T7j, with some specificities 
of its own in the cavity QED version as compared to resonance fluorescence due to the 
incoherent nature of the excitation [18]. Therefore, in this discussion we analyse the 
emission directly from the two-level system; this is the mode to which the sensors are 
now connected. 

A Mollow triplet arises when a two-level system is driven by intense laser light. The 
splitting from very high rungs of the Jaynes-Cummings ladder becomes homogeneous 
giving rise to similar dressed states, |+) = c|0) — s\l) and |— ) = s\0) + c|l), 
with |0) the ground state and |1) the excited state of the two-level system. The 



Two-photon spectra of quantum emitters 

Coherent excitation Incoherent excitation 



17 




Figure 6. Normalised IPS (first row) and 2PS (second and third rows) for the Mollow 
triplet under coherent (left) and incoherent excitation (right). Parameters are chosen 
so that the sidebands are at the same positions Q± = ±6.16g and have the same 
broadening with P a = 2g and (n a ) = 9.5. Vertical/horizontal gridlines mark the 
positions of the upper (green), lower (orange) and central black (peaks), arising from 
the transitions between dressed states depicted in (c). In (d) and (e) the detector 
width, T = 3Pcr/2, is wide enough to filter full spectral peaks. Clear correlations 
appear for the pairs of transitions in (f). In (g) and (h) the detector width, T = P a /2, 
is smaller than the spectral peaks. In this case, the dominant feature is the leapfrog 
triplet of antidiagonal lines, corresponding the to two-photon de-excitation processes 
sketched in (i). 



ladder operator can be decomposed into the four possible transitions, a = |0) (1| = 
° 2 1+) ( — I ~~ cs 1+) (+1 + cs \~) (~ I ~~ s2 \~) (+1- The first and last terms give rise 
to sidebands in the IPS at positions Q± « ±2g^/ (n a ) (at resonance), as showed 
schematically in Fig. [6](c). 

A standard way to describe the laser excitation is through the semiclassical 
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Hamiltonian term H L = Q L (a + a^) that couples the classical term Q L (describing 
a macroscopic laser field) and the quantum term (a + a^) (describing the quantum 
emitter) [T7j. In this approximation, the laser has an infinite coherence time (7L = 0) 
with intensity (n a ) = ^/g 2 and Q± = ±2Q L . This is a successful description of the 
experimental observations, that report the theoretical lineshape, shown in Fig. [6](a). A 
fully quantised picture for the one-emitter laser, such as the Jaynes-Cummings model 
under incoherent pumping, provides the counterpart lineshape in Fig.[6](b). Parameters 
were chosen so that (n a ) coincide in both cases (to 9.5) and the three peaks match 
in position and broadening, with P a = 2g for the Jaynes-Cummings. They are not 
expected to match exactly since in the former case the coherence of the laser is perfect 
and independent of the emitter, while in the latter case, the coherence of the light- 
field is established by the two-level emitter itself. They do share, however, essentially 
the same phenomenology, with the fundamental difference of quantisation or not of the 
light-field between the two. In the lineshapes, this causes the important difference of 
a finite linewidth 7 L for the Rayleigh scattering peak on top of the two-level system 
triplet in the Jaynes-Cummings model, while it is a Dirac S function in the conventional 
Mollow triplet theory due to the infinite coherence time of the laser. We will see that, 
although this is a rather innocent departure in the IPS, it has dramatic consequences 
for the photon statistics. 

We now describe the 2PS of these two cases. Correlations between full, well 
separated peaks of the semiclassical Mollow triplet have been the subject of intense 
experimental [4S1 \M EU and theoretical P21 U31 1501 ISH 1521 031 H51 1531 dEJ efforts (note 
that the Mollow triplet of a cavity QED system in strong-coupling under incoherent 
pumping has still not be reported experimentally at the time of writing). Recently, 
the full 2PS (more specifically a closely related generalised Mandel Q parameter) of 
resonance fluorescence was obtained [16] within the generating function formalism, via 
single-molecule photon counting statistics with spectral resolution. The 2PS for both 
cases are compared in Fig.[HJ They are smoother than in the quantum regime, Fig.[2j and 
with patterns more dominantly aligned along diagonals and antidiagonals rather than 
horizontal and vertical lines, implying that the system is close to a macroscopic, classical 
state as opposed to being dominated by correlations between well defined real states. 
One can understand the main features of these plots again through a decomposition 
into two-photon dressed state transitions: 

™ = s 3 c\-) ( (+| +)-(-| -)) (+| + cV |+) ( (+| +) - (-1 -)) (+| 

- c 3 s |+) ( (+| +) - (-1 -)) (-1 - cV I-) ( (+| +) - (-1 -)) (-1 . (18) 

Each of the 2 3 terms represents a possible deexcitation path and all together, they 
interfere destructively to give rise to the expected total antibunching a a = 0. The first 
line corresponds to the four cascade transitions depicted in Fig. |6](f), that start from 

l+>: 

• (f2 + ,0): One would expect this combination to produce bunched or cascaded 
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Figure 7. Two-photon correlators for frequencies of interest of the Mollow triplets as 
a function of the detector linewidth, T. (a) Correlations of a two-level system driven by 
a large coherent pump (Ol = 1507^) for the four most relevant (cji,^) configurations: 
(0,0) in a solid red line, (f2 + ,f2_) in dashed orange, (f2 + ,f2+) in dotted-dashed blue 
and (fi_|_,0) in dotted blue. Four regions of detection can be defined (separated by 
vertical gridlines). Correlations from the central peak (solid red) exhibit a bunching 
region of [<t](0, 0) ~ 3 at small T when the coherent scattering is a significant 
fraction of the filtered light. Vanishing T yields the statistics of the classical laser 
since the semiclassical zero linewidth cannot display indinstinguishability bunching, 
(b) Comparison of the central peak correlations, (0,0) (solid red), in the case of an 
incoherently pumped two-level system in free space (dotted black) and in a cavity 
(solid black). The laser intensity here is lower: Q± = ±6.16g as in Fig. [6] The linear 
regime, <^ 7<r> is plotted with a dotted red line as a reference in units of j a = g. 

(2) 

When the photon field is quantized, the correct limit J [a] (0,0) = 2 when r — » is 
recovered. 



photons but, as it corresponds to two different paths with amplitudes of probability 
with opposite signs (first two terms in Eq. (18)), destructive interference leads to 
<7p 0) = 0. Such debunching effect occurs within the detector timescale 

i/r [491 EH- 

• Their cascade configuration produces a strong bunching, with a 
well defined time order that depends on the detuning with the laser [T3] . 
However, at resonance, destructive interference debunches again the statistics to 

fl£ 2) [<7](n + ,n_) = i. 

• (0, 0): A similar situation of bunching and interference leads to #r^[0"](O, 0) = 1. 

• (fi +J fi_|_): photons from the same sideband do not form a cascade so their 
simultaneous emission is suppressed (despite the indistinguishability bunching): 

g ^[a](n + ,n + ) = o. 

These ideal correlations are shown in Fig. [7|(a) in the full peak detection region (in 
yellow), where £7l 3> IV With a less intense laser such as than in Fig. [6j the features 
are still visible but slightly tempered due to the overlap between peaks, as is the case 
in experiments [TT] . 

Thanks to the ease of use of our general solution, we are able to compare the 
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semiclassical model for the laser |6j(d) with the Jaynes-Cummings model in the lasing 
regime [6](e). In these qualitatively similar plots, the interference effect described above is 
even more evident as it extends beyond the four points (fl±, 0), forming two blue rings 
of antibunched correlations. These interference rings are more clear in the strongly 
driven situation depicted in Fig. [8j Their origin, as for the leapfrogs, stems from the 
uncertainty introduced by including the detector physics, which extends the interference 
appearing at the points (f2±,0) and (fi±,£2±) to two circles. At the center of these 
interference rings, (Q±/2, f2±/2), two red spots of enhanced emission appear from the 
concentration of leapfrog outer lines in the Jaynes-Cummings model, Eqs. (22). With a 
narrower detector, leapfrogs are enhanced due to the longer uncertainty in the time of 
two-photon emission, as shown in Fig. [6](g) and (h). In the same way that the Jaynes- 
Cummings multiplets from high rungs of excitation converge to the Mollow triplet in 
the IPS, the antidiagonal leapfrog lines in the 2PS converge to a leapfrog triplet at: 

cji + uj 2 ~ , cji + uj 2 ~ Q± , (19) 

as schematically depicted in Fig. |6](i). They are more pronounced in Fig. |6](g), as the 
perfect laser approximation, Hl, generates dressed states homogeneously split for all 
intensities f^- 

If the detectors are narrower than the peaks, r < T ai the interference effect 
described above disappears (see Fig.[7|. Furthermore, as has been previously discussed, 
when T is smaller than any peak width, photons are detected from any point of 
the dynamics and appear as uncorrelated. For different frequencies this means 
limr^o M (^i ' ^2) = 1 while for equal frequencies, limr^o 9r\ (J ]( UJ i UJ ) = 2 due to 
the indistinguishability of the photons. In this limit, the statistics observed can no 
longer be attributed to the system, but to the detection process. The two-level system 
excited incoherently (with or without the cavity) always recovers this limit correctly as, 
in the Jaynes-Cummings model, the smallest width, 7 L is still finite and sets a clear lower 
boundary (see Fig. [71(b)). In contrast, we have a totally different result under coherent 
excitation of a perfect laser whenever uj\ + uj 2 = 2cjl = 0: limr-^o 9r 0) = 1 and 

limr^o 9r\ a ](^+^ ^-) = 2- Due to the implicit assumption of zero-linewidth for the 
laser field, it is not possible in these cases for the detector line width to be thinner and 
reach the uncorrelated limit. Instead, the limit r isolates completely the Rayleigh 
peak so correlations become exactly 1, i.e., those of the laser. This occurs with a detector 
linewidth below r min , defined as that at which the coherent part of the filtered IPS is 
as large as the incoherent part, 

CM(o)-C h M(o) 5 (2°) 

where S t [<t](u) = S^ oh [a](uj) + S^ coh [a](uj) . In the limit of intense lasing, r min « 
1^/(4^). In the intermediate region, r min < T < T a: correlations from the central 
peak are bunched (and equal to 3 in the ideal case), due to the mixture of scattered and 
emitted light. 

In real lasers, 7l 7^ due to phase fluctuations. Non-monochromatic theoretical 
models, where the phase varies stochastically [55j [56], take into account both the 
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amplitude and speed of fluctuations, recovering not only a finite 7l for the laser 
mode [53, [55] but also physical correlations at all limits [58J. The Jaynes-Cummings 
model or one-atom laser is free from this pathology as fluctuations are intrinsic to the 
dynamics, provided by the interplay between coherent exchange, incoherent pump and 
decay. However, not being fully coherent unless very far above threshold, one requires 
a very good system (j a <C g) in a very intense and well defined lasing regime, where 
7l < 7a, r min , that is P a > y/2g. 

Like for any other system, very broad filters T 3> ^+ recover the full antibunching of 
the two-level system, therefore limr-K»#r 2 ^[°"](0j 0) = fl^M = 0- This is independent of 
the type of excitation. Under very low coherent excitation (dotted red line in Fig. [7|(b)), 
where the IPS is dominated by the coherent part, correlations monotonously go from 1 
to as 




This regime has been exploited to create an ultra-coherent single photon source, i.e., 
indistinguishable in frequency thanks to the inherited long laser coherence time [59l l60l 
I6T1 1621 163] . It is interesting to note that filtering the coherent peak, T < 7^, yields 
the coherent statistics of a laser and to the destruction of the antibunching even if the 
coherent part strongly dominates the IPS. Only the full emission has the property of 
interfering destructively and provide photons one by one. 

Finally, the impact of detuning from the laser (in the semiclassical approximation) 
is shown in Fig. [8} In the absence of pure dephasing, the IPS remains symmetric around 
the laser frequency (here uj^ = 0). Interestingly, the 2PS does not, even if detuning is 
small as compared to the driving, as shown in (b). The asymmetry in the 2PS consists 
in the disappearance of the ring of interferences nearer to the two-level system frequency 
(in this case uj g > uj^). When detuning dominates and Q+ w u) a , as in (c), the region 
around this peak becomes more antibunched, as it correspond to a real transition in a 
two-level system. On the other hand, the region around the other sideband, (fl_,f2_), 
becomes more bunched, as it corresponds to a virtual laser transition [12l [13]. Two 
red lines cross at only interrupted by the blue interference ring. A strong 

bunching point emerges at (f2 +5 fi_), accompanied by the typical cascade r-dynamics, 
making the system suitable for heralded single photon emission [TT] . 

One can understand the asymmetry in the 2PS with detuning, g^[a](Q + , Q+) ^ 
g^\o]{Vt_ ) fi_), as a manifestation of the dissimilar two-photon dynamics of the dressed 
states, namely, the r-dynamics of the colour-blind correlations corresponding to each 
transition: g^[\ — ) (+|](t) and fl^[|+) (— |](t). These two functions are identical in 
the limiting r- values, equal to at r = and to 1 at r — » 00, but evolve differently 
at intermediate r. In the limit of high excitation and full filtering of each peak, such 
dressed state two-photon dynamics is overall well mapped by the frequency-resolved 
correlations, g^[a](Q +) Q + , r) and ^p 2 ^[a](^_, r), as long as r > 1/r. However, as 
we have shown, the r = value is not well mapped, being larger than zero when the 
detection is included (spoiling perfect antibunching). The frequency-resolved functions 
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Figure 8. 2PS for a two-level system under coherent excitation deep in the Mollow 
triplet regime (0 + = y/ (2^l) 2 + — ^l) 2 = 300)^) in the cases where (a) the laser 
is on resonance with the emitter (Ol = 1507^, as in Fig. (7^a)) and (b) the laser is 
slightly detuned (Ol = 1007^) or (c) largely detuned ($\, = 107 CT ). In contrast with 
the IPS, the 2PS is asymmetric with detuning. Also, in this regime where a highly 
coherent field is involved, a new type of pattern emerges in the correlations, namely, 
circles. 



at r = 0, include an average over the short-time dressed-state dynamics, within the 
detector time-scale, 1/r. Therefore, the interaction with the filter/detector can be 
considered a dephasing mechanism of the dressed-state dynamics, which breaks the 
symmetry of the 2PS even at r = 0. This view is confirmed by the behaviour of 
the physical time-dependent IPS [21], in the absence of pure dephasing, which is also 
asymmetric in general before reaching its symmetric steady state lineshape (not shown). 

6. Conclusions and further directions 

In conclusion, we have performed the first systematic investigation of the two-photon 
spectra (2PS) of a variety of quantum optical systems, using our recently developed 
formalism to compute frequency and time resolved TV-photon correlations pQ. We have 
focused here on the case N = 2 of two-photon correlations at zero time delay. Thanks to 
our formalism, such investigations can be generalised to higher photon number and/or 
arbitrary time delays. We studied systems of increasing complexity, starting from the 
simplest possible case that is the quantum harmonic oscillator and proceeding with 
the two-level system and then all their possible combinations, namely, two coupled 
harmonic oscillators, two two-level systems and a mixture of both that amounts to the 
celebrated Jaynes-Cummings model. The latter, in contrast with our starting point, is 
an extremely rich and complicated system that stands as the pillar of cavity quantum 
electrodynamics. We have outlined what are the common features shared by all these 
systems and what are those specific to their underlying structure and dynamics. These 
results constitute the backbone for two-photon spectroscopy. The main finding can be 
summarised as follows: 
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There is a universal bunching of photons when filtering an emission line below 
its line width, regardless of its inherent statistics, due to the indistinguishable and 
uncorrelated arrival of photons (indistinguishability bunching). This manifests as a 
diagonal on all 2PS. 

An antibunching arises for a two-level system with a characteristic butterfly-shape 
that is reproduced by any two-level transition that can be resolved in isolation 
from a more complex level structure (see Ref. [64] for an example within a four- 
level system). 

The 2PS of various types of coupled systems differ greatly from each other even 
when their IPS are identical, and within the same system, some symmetries of 
the IPS are lifted in 2PS. This shows the much higher degree of characterization 
accessible via two-photon spectroscopy. 

The 2PS of coupled harmonic oscillators features an anticrossing (two hyperbolas) 
of uncorrelated emission corresponding to the independent and classical polaritons 
dynamics. Coupled two-level systems or the Jaynes-Cummings model on the 
other hand give rise to rich patterns of correlations permitting to investigate the 
underlying physical picture in terms of relaxation between quantum states. 

Such correlations are due to fundamental processes that can be identified through 
simple equations that locate them in the 2PS. They are i) leapfrog processes, i.e., 
two-photon emission jumping over an intermediate rung of the level structure 
through a virtual state of indeterminate frequency and ii) polariton-to-virtual- 
state ant icorrelat ions, where the emission does not correspond to dressed states 
transitions but take into account the dynamical nature of the system, with the 
final state provided by a virtual state or fluctuation of the system. 

Leapfrog processes are particularly noteworthy for applications in quantum 
information processing. They are present in any nonlinear quantum system, as 
simple as an anharmonic oscillator. TV-photon leapfrogs with N > 2 also take 
place, though the higher the TV, the weaker the process. In the Jaynes-Cummings 
system, they are evidenced by higher order correlations, g^fai, • • • ,wn), a t the 
frequencies 

N 

Y,u k = Nuj a ±(E n -E n _ N ) or (22) 
i=i 

= Nu a ± (E n + E n _ N ) , (23) 
but their further characterization is out of scope of this text. 

The intricate pattern of TV-photon correlations that is formed by the combination of 
all the above processes in the Jaynes-Cummings model, that is neatly resolved in a 
system well into the strong-coupling regime, evolves when brought in the lasing 
regime into the qualitatively different 2PS of the Mollow triplet, that features 
new types of patterns such as circles rather than simply straight lines. For 
high rungs of excitation, n 3> 1, the lines in Eq. (22) converge to the central 
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antidiagonal, J^c^ « 0, while the lines in Eq. (23) agglomerate at the outer 
positions, J^o^ ~ ±2y/ng. This formation of the Mollow triplet in the IPS gives 
rise to leapfrog triplet in the 2PS, with the same origin as in the linear regime in 
terms of two-photon transitions between dressed states. 

Beyond spelling out the dynamics of emission at the quantum level, the theory of 
two-photon spectroscopy also allows to address fundamental theoretical issues of the 
quantum formalism and illustrates the deep link between the quantum dynamics and 
the detection process. For instance, in the limit of an ideal detector, the semiclassical 
description of the Mollow triplet fails to recover the fundamental indistinguishability 
bunching. This is due to the artificial 5 line of the scattering peak produced by the 
semiclassical approximation. This shortcoming can be solved by turning to a fully 
quantized theory or upgrading the semiclassical theory to get rid of the artifacts caused 
by a fixed phase. 

Another conclusion from our study is that the dynamics of "real processes", 
involving photon by photon de-excitations, are related to the system parameters and 
less influenced by detection. The corresponding correlations are best resolved when the 
related peaks are separated and fully filtered. On the other hand, "virtual processes" , 
involving virtual states such as the leapfrogs and polariton-to- virtual state interferences, 
happen within the time of interaction with the detector, 1/r, and conserving energy 
within its linewidth, T. Therefore, their correlations become more prominent using 
narrow filters. 

We have thus amply demonstrated that two-photon spectroscopy unravels a rich 
two-photon dynamics and interference effects in open quantum systems by a precise 
disposition of filters of given resolutions. Such correlations can be further taken 
advantage of by considering finite time delays and/or optimising the frequency windows. 
Further interesting systems to apply two-photon spectroscopy are ultra strong coupling 
systems [65], closely spaced atoms in optical lattice [66], the biexciton two-photon 
emission in a quantum dot f42l 164] or the dynamics of Bose-Einstein condensation [67J . 
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